Anomalous near-field heat transfer between a cylinder and a perforated surface 
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We predict that the radiative heat-transfer rate between a cylinder and a perforated surface depends 
wtwmonotonically on their separation. This anomalous behavior, which arises due to near-field effects, is ex- 
plained using a heuristic model based on the interaction of a dipole with a plate. We show that nonmonotonicity 
depends not only on geometry and temperature but also on material dispersion — for micron and submicron 
objects, nonmonotonicity is present in polar dielectrics but absent in metals with small skin depths. 



PACS numbers: 

Disconnected bodies of different temperatures can ex- 
change energy through stochastic electromagnetic waves [1- 
7], a phenomenon known as radiative heat transfer that un- 
derlies many naturally occurring and technologically rele- 
vant processes [7] . Recent advances in microfabrication and 
metrology have enabled experiments that can now regularly 
probe this phenomenon at micron and submicron scales [8, 9]. 
At such small scales, unusual near-field interactions arise [7], 
but non-planar geometries in this context are only just begin- 
ning to be explored [10-17, 19, 20]. In this Letter, inspired 
by our previous work on Casimir repulsion [18], we demon- 
strate that the heat transfer between a cylinder (or elongated 
object) and a perforated surface (e.g. a ring) can vary non- 
monotonic ally with respect to their mutual separation, in con- 
trast to what has been observed in all previous geometries [7] . 
This anomalous effect stems primarily from the contribution 
of dipolar (near) fields, as we illustrate by a heuristic model in 
which the cylinder is modeled as a quasi- static dipole and the 
ring as an infinitesimally thin plate with a hole. We find that 
nonmonotonicity weakens (and eventually disappears) when- 
ever the geometrical and material parameters of the objects 
deviate significantly from the dipolar regime: for cylinders of 
equal or nearly equal aspect ratio, large ring thicknesses, large 
temperatures, or metals with small skin depths (such as gold), 
the usual monotonic dependence is observed. We note that, in 
contrast to conventional geometries, this effect cannot be pre- 
dicted even qualitatively by "additive" approximations such 
as the well-known proximity approximation. 

In the far field (object separations d much greater than the 
thermal wavelength At = hc/ksT), radiative heat transfer is 
dominated by the exchange of propagating waves and is thus 
nearly insensitive to changes in separations. In the (less stud- 
ied) near field (object separations d < At), not only are inter- 
ference effects important, but otherwise-negligible evanescent 
waves also contribute flux [6, 7]. Such near-field effects have 
been most commonly studied in planar geometries, where they 
are known to lead to monotonically increasing heat transfer 
rates with decreasing d, resulting in orders-of-magnitude en- 
hancements of the total heat transfer (which can even exceed 
the far-field black-body limit at sub-micron separations [7]). 



Thus far, little is known about the heat transfer characteristics 
of bodies whose shapes differ significantly from the planar, 
unpatterned structures of the past. Recent theoretical progress 
include predictions for a handful of new geometries, includ- 
ing spheres [10, 14] and cones [19] suspended above slabs, as 
well as patterned surfaces [11-13, 15, 20]. Our work extends 
these studies to yet another class of possible geometries: in- 
terleaved bodies whose near-field interactions and shapes lead 
to anomalous heat-transport phenomena. 

The heat transfer rate H between two objects held at tem- 
peratures Ti and T 2 can be expressed in the form [6, 7]: 

rOO 

H= duj [G(u, Ti) - G(u, T 2 )} $(cj), (1) 
Jo 

where is the flux spectrum (the time-averaged flux into ob- 
ject 2 due to current sources in object 1), and 6(cj,T) = 
hu/[ex.p(hbj/kBT) — 1] is the mean Planck energy per oscil- 
lator at frequency uj and temperature T. (Note that = 1 for 
black bodies that capture all of one another's radiation.) We 
compute by exploiting two recent computational methods: 
a fluctuating surface-current (FSC) formulation involving the 
solution of an integral equation at each frequency [17], and a 
Langevin finite-difference time-domain (FDTD) formulation 
in which one explicitly time-evolves Maxwell's equations in 
response to (broad-bandwidth) stochastic sources inside the 
bodies [12]. In order to distinguish the effects of geome- 
try from those of material dispersion, we begin by consider- 
ing a simple model material: a lossy dielectric with a broad 
(low-dispersion) absorption peak, given by e(uj) = — 
ct/(uJq — uj 2 — ijuu), with = 12.5, a = 4 x 10 2 (c//im) 2 , 
ujo = 0, and 7 = 60 (c//im), corresponding to roughly- 
constant Ree « 12 and large Ime > 1 over relevant fre- 
quencies. Later, we consider realistic materials and show that 
material dispersion also plays a crucial role. 

We previously studied the radiation of isolated cylinders 
and rings [17]. In this work, we consider the new phenom- 
ena that arise when these two objects are brought into close 
proximity so that near-field effects are present. The dashed 
green lines in the inset of Fig. 1 show the flux spectrum of 
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FIG. 1: Heat-transfer rate H from a room temperature ring to a 
cylinder (or sphere) of fixed radius R = 0.1/xm and aspect ratio 
A = L/2R held at T — 0, as a function of their center-center sep- 
aration d. H is normalized by the heat radiation of the isolated ring 
H r . Both objects are lossy dielectrics with Res « 12 (see text). In- 
set shows the flux spectrum <3>(u;) of the A = 5 configuration at three 
separations, and also of the isolated cylinder and ring (dashed lines). 



an isolated cylinder (<1> C ) of radius R = 0.1/im and aspect 
ratio A = L/2R = 5, and of an isolated ring (<l> r ) of outer 
(inner) diameter D = 2 /mi (W = 0.8/im) and thickness 
h = 0.05D, as a function of frequency (units of 2R/X). In 
this long-wavelength limit (A ^> R), large Ime and the ab- 
sence of geometric and material resonances means that both 
objects emit significantly less than an ideal black body, so that 
$ c < <$>r < 1. 

The inset in Fig. 1 also shows the flux spectrum when 
the two objects are brought into close proximity (solid lines), 
at three different center-center separations d, showing dra- 
matic changes from the isolated case. is computed via 
both FDTD (noisy curves) and FSC (smooth curves) meth- 
ods, showing excellent agreement; the remaining calculations 
use FSC only. Compared to isolated objects, the increasing 
contribution of evanescent waves in both objects leads to an 
overall increase in the flux at low frequencies. [Note that the 
peak in $ at A « 10 _3 i? is a consequence of material disper- 
sion: the loss tangent of the material ~ Ime/ Re e —> oo as 
A — » oo, leading to zero radiation; the peak in the spectrum 
occurs at the cross-over wavelength for which Ime ~ Ree. 
Low-frequency cut-offs in the near-field enhancement occur 
in highly conductive materials, such as gold (below).] Most 
interestingly, however, the enhancement in <1> here does not 
increase monotonically with decreasing d: increases from 
d = oo to d « 0.4L, but decreases from d « 0.4L to d = 0. 

To explore the geometry dependence of this near-field be- 



havior, we now examine the overall heat-transfer rate H as 
a function of d instead of the spectrum. In particular, Fig. 1 
shows H from a room-temperature ring to a cylinder at T = 0, 
for multiple aspect ratios A. For comparison, H is normal- 
ized to the radiation rate of the isolated ring H Y . For large 
anisotropy A = 5, H first increases as the two objects ap- 
proach each other due to the usual near-field enhancement, 
and then decreases as d — > 0, peaking at a critical separa- 
tion d c « 0.4L. Unlike previously studied structures involv- 
ing non-interleaved objects, the heat transfer in this geome- 
try does not diverge as d — » 0: although the two objects 
approach each other in this limit, they never touch. Also, 
H — » as d — » oo due to the finite size of the two objects. 
As A decreases (keeping R fixed), corresponding to increas- 
ingly isotropic cylinders, the nonmonotonicity becomes less 
pronounced, and is completely absent in both the A = 1 con- 
figuration (small anisotropy) and for a sphere (open red cir- 
cles). As expected, there is an overall decrease in H with de- 
creasing L due to the decreasing volume of the cylinder. Non- 
monotonicity also slowly disappears as the ring thickness h in- 
creases to h ~ 0.5L, leaving a relatively wide range of thick- 
nesses over which the effect can be observed. Furthermore, as 
expected, the strength of the nonmonotonicity H(d c )/H(0) 
grows larger as the cylinder surface approaches the rim of the 
ring (corresponding to larger R or smaller W) due to near- 
field effects, and also for larger D due to the larger surface 
area. More interesting however, is the fact that nonmonotonic- 
ity persists even when the cylinders are shifted laterally (shifts 
< 0.5VF), an asymmetric configuration that is likely to occur 
in experiments. Finally, we find that H(d c )/H(0) increases 
as h, i?, L — » (for fixed A), as seen below (Fig. 2). 

Cylindrical symmetry allows us to decompose into az- 
imuthal angular components m (fields ~ e zm6> ), implemented 
in FDTD with cylindrical coordinates. Our calculations reveal 
(not shown) that most (though not all) of the nonmonotonic 
dependence comes from the contribution of dipolar (m = 0) 
fields, which dominate the heat transfer at these long thermal 
wavelengths At ^ R,L. At these wavelengths, a cylinder 
with large A will act like a fluctuating dipole oriented mainly 
along the symmetry axis. Since the fields generated by a fluc- 
tuating dipole are polarized mostly along the dipole axis and 
since current fluctuations in the thin ring are polarized mostly 
along the plane of the ring, it follows that the fields induced 
by a long cylinder will do less work on currents in the ring 
whenever the objects are nearly co-planar (d —> 0). 

We quantify this argument by focusing on a simple (al- 
beit heuristic) model in which the cylinder is modeled as a 
dipole of electric polarizability cx E (u) and the ring as an 
infinitesimally thin plate with a hole. For convenience, we 
only consider the non-retarded (quasi-static) limit of small 
dipole separations d <C At, in which case the heat trans- 
fer rate between the dipole and plate can be expressed as 
H d = jduj [6(o;,Ti) - e(o;,T 2 )] $dM> where the flux 
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FIG. 2: Heat-transfer rate if between the ring and cylinder of Fig. 1, 
for fixed cylinder aspect ratio A = 5 and thin ring thickness h — 
5 x 10 _3 D, as a function of d and for multiple cylinder lengths L. 
H is normalized by the radiation of the isolated ring H r multiplied 
by the cylinder area A c . The shaded region denotes separations over 
which the two objects are interleaved. Top inset shows the critical 
separation d c of largest if as a function of L, with shaded areas again 
denoting interleaved configurations. Bottom inset shows H in the 
limit h,R, L —> 0, as computed by a heuristic model (see text). 



spectrum is given by [5, 21, 22]: 



2 ^ 

$ d (w) = - Vw 2 Im[af(w)] Im[G«(r d ,r d ) 

7T f * 



(2) 



Here, Gij (r, r') is the electric Green's function of the plate — 
the electric field in the zth direction at r due to a dipole source 
in the jth direction at r' — and rd is the location of the dipole. 
The calculation of G^ for a perfectly conducting, infinitely 
thin plate was carried out in Ref . 23 for the purpose of comput- 
ing the Casimir-Polder force in that idealized system. How- 
ever, because perect conductors do not radiate (Im G^ = 0), 
we merely exploit that expression as the starting point of a 
quasi-static perturbative calculation in which the plate is as- 
sumed to have a small amount of absorption. In particular, 
we are interested in computing the dissipated power or Ohmic 
losses on a plate with small but finite conductivity a, given the 
quasi-static fields at the surface of the perfectly conducting 
plate. Following Ref. 24, the resistive losses on the plate are 
~ ImG^(r d ,r d ) ~ / d 2 v a(r)|G°(r, r d )| 2 , where G% is 
the electric Dyadic Green's function of the unperturbed (per- 
fectly conducting) plate, and the integral is performed over 
the plate surface. It turns out that G° exhibits non-integrable 
singularities at the rim of the hole, a well-known artifact of 
the idealized nature of corners and wedges in electromag- 



netism [24]. While the electromagnetic energy corresponding 
to these fields is finite when integrated over all space, the sin- 
gularity in the fields is problematic for the perturbation theory 
since the form of the perturbation considered here requires 
that G° be integrated only over the plane of the plate. Essen- 
tially, this model does not account for the finite thickness of 
the plate, and consequently fails to capture effects associated 
with the finite penetration or skin-depth S = c/(lj Im y/e) of 
fields. Therefore, we use S as a cutoff lengthscale to regularize 
the integral near the rim of the hole. Surprisingly, and despite 
its many shortcomings, this heuristic model captures most of 
the features of interest. 

The bottom inset in Fig. 2 shows this heuristic model's heat 
transfer rate if d (in arbitrary units) from a room-temperature 
plate of conductivity a = ujlms, infinitesimally small thick- 
ness h — » 0, and hole diameter W, to a small cylinder of 
radius R <^W and electric polarizability ol e held at T = 0, 
as a function of their separation d. is computed by Eq. (2) 
using the dipole model above, with the polarizability of the 
dipole taken to be that of a uniform prolate spheroid [22] of 
aspect ratio A = L/2R and permittivity e (same as above). 
For A = 5, corresponding to a highly anisotropic object, one 
observes the expected nonmonotonic behavior, with the criti- 
cal separation d c « 0.3W now determined by W. As before, 
nonmonotonicity decreases with decreasing A, disappearing 
completely in the limit A — » 1 of an isotropic (spherical) ob- 
ject with polarizability a E = 4/37i\R 3 (£ — l)/(e+2). In com- 
parison with Fig. 1, we observe that the simple model quan- 
titatively captures the onset of nonmonotonicity at A > 2. 
Since the model represents a point-dipole limit, we also com- 
pare to exact calculations for fixed A = 5 by letting L — > 
and h = bxlO~ 3 D, as shown in Fig. 2. While nonmonotonic- 
ity is present for all L, the scaling of the critical separation d c 
with L changes qualitatively as L — » 0. Specifically, as shown 
by the top inset of Fig. 2, d c « 0.3W for L <C W, in quanti- 
tative agreement with the dipole model, while d c ~ 0.4L for 
L ^> W, with the cross-over regime occurring for L w 0.5 W. 
We find that for L < 0.SW, the onset of nonmonotonicity (d c ) 
occurs before the two objects are interleaved, i.e. when there 
is a separating plane between the objects. 

Aside from geometry, changes coming from either temper- 
ature or material dispersion leading to deviations from the 
ideal dipole regime can also weaken nonmonotonic behavior. 
Thus far, we have restricted ourselves to studying micro- scale 
bodies near room temperature (which emit preferentially at 
infrared frequencies), corresponding to large thermal wave- 
lengths At ^ R. At larger temperatures (At < R) however, 
such bodies can no longer be well described as dipole emit- 
ters, and thus <£> no longer exhibits nonmonotonic behavior. 
We find that nonmonotonicity persists at temperatures well 
beyond the mere T « 300 K considered here, so long as the 
feature sizes of the objects involved remain at or below the 
micron scale. Realistic materials often exhibit substantial ma- 
terial dispersion at or near infrared wavelengths, and this can 
also significantly alter the dipole picture above. This situa- 
tion is depicted in Fig. 3, which shows H for various ma- 
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FIG. 3: Heat transfer rate H between the ring and cylinder of Fig. 1, 
for fixed cylinder aspect ratio A = 5, as a function of separation d. H 
is normalized to the isolated cylinder radiation rate H T , and plotted 
for multiple material configurations. Insets show the corresponding 
flux spectra $(cj) at three different (labeled) separations. 



terial configurations, including metals (gold and indium tin 
oxide) and polar dielectrics (doped silicon). (The Au and 
ITO dispersions are determined by Drude models with plasma 
frequencies uj p = 1.367 x 10 16 rad/s and uj p = 1.4739 x 
10 15 rad/s, and relaxation rates 7 = 5.317 x 10 13 rad/s and 
7 = 1.5347 x 10 14 rad/s, respectively, whereas the doped sili- 
con dispersion is given by the model of Ref. 25.) Interestingly, 
we find that nonmonotonicity is completely absent for the 
Au and highly-doped silicon (10 22 cm -3 ) configurations, and 
present in the ITO and lesser-doped silicon (5 x 10 18 cm -3 ) 
configurations. The reason for the discrepancy comes from 
the fact that highly conductive metals and polar dielectrics re- 
spond very differently to incident light. In particular, for met- 
als, the electric dipole approximation (above) breaks down 
for skin depths S <C R, L: in that limit, eddy currents in- 
duced on the surface of the metallic objects also lead to large 
magnetic dipole moments [21, 22], with magnetic polarizabil- 
ities ol h ^> ol e . Unfortunately, the interaction between a 
magnetic dipole moment and a plate does not exhibit the de- 
sired nonmonotonic effect (at least in this geometry), which 
explains the results in the case of Au and highly-doped silicon 
cylinders, whose skin depths 5 « 10 _2 /im <C R at infrared 
wavelengths. If we scale the entire structure down to much 
smaller scales R,L <C S (not shown), we find that nonmono- 
tonicity is restored. It follows that for highly conductive ma- 
terials, one obtains the desired nonmonotonic effect only for 
h,R,L<^S < A T . 

The insets in Fig. 3 show the flux spectra of the various ma- 
terial configurations at three separations (labeled A, B, and 



C). The flux spectrum for Au and highly-doped silicon are 
monotonic with d at all frequencies. On the other hand, ITO 
and lesser-doped silicon are nonmonotonic with d only at cer- 
tain frequencies. For any material, is monotonic at high 
frequencies (A <C R, L), where the cylinder no longer acts 
as a dipole, but this transition occurs even more rapidly for 
ITO (which is only nonmonotonic for low frequencies) be- 
cause of the aforementioned skin-depth effect. For silicon, the 
situation is greatly complicated by the presence of multiple 
geometric resonances (e has a single absorption peak), aris- 
ing from the ability of the fields to probe the interior as well 
as the surface of the dielectric. It turns out that only two of 
the resonant peaks in this case exhibit nonmonotonicity with 
d, corresponding to resonances with the necessary dipole-like 
polarizations. 

Similar and even more pronounced nonmonotonic behav- 
iors should arise in other geometries, so long as the suspended 
objects (regardless of shape) are sufficiently anisotropic (be- 
have dipole-like) and radiate primarily in the direction orthog- 
onal to the patterned surface. An interesting structure to ex- 
plore in the future is a nanowire array suspended above a pe- 
riodically patterned thin film. 
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